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PHYSICS FUNDAMENTALS OF LUMINOUS 
ACCRETION DISKS AROUND BLACK HOLES 


O. M. Blaes 


Abstract 

These lectures provide an overview of the theory of accretion disks 
with application to bright sources containing black holes. I focus on 
the fundamental physics of these flows, stressing modern develop¬ 
ments and outstanding questions wherever possible. After a review 
of standard Shakura-Sunyaev based models and their problems and 
uncertainties, I describe the basic principles that determine the over¬ 
all spectral energy distribution produced by the flow. I then describe 
the physics of angular momentum transport in black hole accretion 
disks, stressing the important role of magnetic fields. Finally, I discuss 
the physics of radiation magnetohydrodynamics and how it might af¬ 
fect the overall flow structure in the innermost regions near the black 
hole. 


1 Introduction 


Accretion disk theory was developed in a flurry of activity in the early 
1970’s in a series of important, seminal papers |95 , 102| , 88|, |75|, [72|, p0||38|, 
|103| |. It has met with considerable success in explaining the observations 
we see in a number of classes of accretion-powered sources, particularly 
cataclysmic variables. However, it does a rather poor job of predicting and 
explaining the observations of black hole sources, both in the context of X- 
ray binaries and active galactic nuclei. This is particularly disappointing in 
view of the fact that much of the early development of the theory was done 
with this application in mind. Accretion disk theory has always suffered 
from a severe flaw: angular momentum transport and energy dissipation 
in observed flows must be due to nonlinear physics (“turbulence”), and 
this problem is swept under the carpet by parameterizing it away in terms 
of an anomalous “viscosity”. Over the past thirty years, this sad state 
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of affairs has persisted, although a major breakthrough was made in 1991 
with the discovery that the magnetorotational instability (MRI) might be a 
generic source of turbulence in disks ||l6|,^2). Thanks to this insight and the 
increasing sophistication and capabilities of numerical simulations, we are 
now beginning to address virtually all aspects of accretion disk theory from 
first principles physics. The coming decade should see such models applied 
to observation, finally getting us to the point where we are able to test the 
basic physics, and not fudge factors. We are therefore living and working 
in an exciting time! 

In these lectures, I will give an overview of the current state of the 
theory. Although I will be discussing accretion disks in some generality, 
I will for the most part concentrate on disks around black holes. Unless 
explicitly noted, I will always use Newtonian physics, even though the flows 
around black holes necessarily involve relativistic effects, both special and 
general. I do this purely for the sake of pedagogy. Relativity is not difficult 
- in fact, it is the easy part of black hole accretion disk physics. As we 
will see, all the really nasty (albeit interesting!) issues lie in the radiation 
and plasma physics, and the basic principles involved here are most quickly 
understood within a Newtonian framework. Almost everything I discuss 
has been worked out in full general relativity somewhere in the literature, 
and I give references to such work wherever possible. 


Many good textbooks exist which discuss aspects of the basic physics 
of accretion disks, e.g. 42,67,1041, and the reader is encouraged to consult 
them for a more elaborate treatment of the now-standard ideas. As I already 
mentioned, however, many of these ideas will soon be superseded by much 
more sophisticated models based on real physics. The best overview of the 
MRI is the review by Balbus & Hawley [|n|, although some advances have 
taken place since that article was written. 


I will begin in section 2 with a brief overview of standard models based 
on the Shakura-Sunyaev a-prescription for the anomalous stress that must 
exist in these flows, as these are still the models that people use for compar¬ 
ison to observational data. I will focus on recent developments in this area, 
including the role of advection in “slim disks”, the role of torques exerted 
across the innermost stable circular orbit, and the vertical transfer of me¬ 
chanical energy out of the disk into a corona. In section 3,1 will discuss the 
calculation of photon spectra from these models. I will then return to first 
principles in section 4 by examining the basic origin of anomalous angu¬ 
lar momentum transport in disks, emphasizing the role of magnetic fields. 
Finally, in section 5, I will discuss the role that radiation magnetohydro¬ 
dynamics may play in determining the dynamics and thermodynamics of 
the innermost regions of the flow, where most of the observed radiation is 
thought to originate. 
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2 Shakura-Sunyaev Based Models 

Most, but not all, models of accretion-powered sources assume that rotation 
is an important source of dynamical support in the accretion flow. Gas ac¬ 
creting onto a compact object from an orbiting companion star is endowed 
with considerable angular momentum arising from the orbital motion of 
the binary itself. The huge dynamic range in radii between the fuel source 
and the radius of the central object in young stellar objects (YSO’s) and 
active galactic nuclei (AGN) also strongly suggests that rotation is impor¬ 
tant. If material at large radii has some nonzero angular momentum with 
respect to the central object, and if that angular momentum is conserved 
on the way in, then the material is likely to become rotationally supported. 
This is simply because centripetal acceleration then depends on distance R 
from the rotation axis as R~ 3 , steeper than the R~ 2 dependence of grav¬ 
itational acceleration. (This statement neglects general relativistic effects. 
Sufficiently close to the event horizon, the gravitational acceleration will 
dominate all rotational support and even material with conserved angular 
momentum will flow inward all the way to the singularity.) Direct obser¬ 
vational support for the existence of rotationally supported accretion flow 
structures exists in many images of YSO disks and in the maser disk of the 
low luminosity AGN NGC 4258 (82). 

Angular momentum will generally not be conserved in the flow, and the 
mechanisms of outward angular momentum transport and how they are 
connected to the conversion of orbital mechanical energy into other forms 
are central to understanding how accretion power works. Ordinary fluid 
viscosity is far too weak to be a significant factor, and something else must 
exert torques (either internal or external, or both) on the flow. 

Virtually all models that actually attempt to make contact with the 
observations are based on a phenomenological prescription of a (vertically- 
averaged) anomalous internal stress introduced by Shakura & Sunyaev in 


1973 (l02j| : 


( 2 . 1 ) 


TR 4 , = aP, 


where P is a vertically-averaged pressure. Many workers in the held con¬ 
tinue to think about this stress as a form of “viscosity”, and the stress 
prescription is often explicitly introduced as exactly that! (As we shall see 
below, such naive thinking is very dangerous and can be extremely mislead¬ 
ing.) For example, a commonly used modification of the Shakura-Sunyaev 
prescription is to assume a kinematic viscosity of the form 



( 2 . 2 ) 


where c 2 — P/p is the vertically averaged sound speed (and p is a vertically 
averaged density), H = H ( R ) is the vertical half-thickness of the disk, and 
Dk is the angular velocity of test-particle circular orbits (the “Keplerian” 
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angular velocity). When inserted into the standard viscous form of the 
stress, this gives 

dfi R dtt 

T «* = p,,R dR =aP nr<dR- (2 - 3) 

Here = Sl(i?) is the actual angular velocity in the flow, which may differ 
from fix if the flow is not completely geometrically thin. If ~ fix, then 
within factors of order unity, equation (2.3) gives the same stress as equation 
& The only difference is that the stress now depends explicitly on the 
shear in the flow, just as an ordinary viscous stress would. This simple 
change completely alters the mathematical character of the critical point 
problem in stationary advective flows |H|], but it is not clear that any of 
this is real. 

Yet another uncertainty that has plagued black hole accretion disk mod¬ 
els with Shakura-Sunyaev stress prescriptions from the very beginning is 
what pressure (gas or radiation or...?) to stick into equations (0) or CD 
when radiation pressure is comparable to or greater than gas pressure. This 
occurs in the central parts of standard disk models arou nd b l ack holes ac¬ 
creting at anywhere near the Eddington rate (see eqs. |5.11-[5.2 below), 
right where most of the power is generated. Moreover, the first thing that 
one would try (just inserting the total, gas plus radiation, pressure) leads 
to a disk tha t ca nnot be stationary because it is thermally and “viscously” 
unstable [j72, |l03 l . 

Let us briefly trace the argument for the radiation pressure driven ther¬ 
mal instability, which acts on a much faster time scale ~ 1/ (af l) in geo¬ 
metrically thin disks than the “viscous” instability which acts on the radial 
flow time scale ~ R 2 /(H 2 aQ). Assuming that cooling of the disk proceeds 
through radiative diffusion, the local emergent flux at radius R is given by 


F~ 


acT 4 


(2.4) 


where a is the radiation density constant, c is the speed of light, T is a 
measure of the disk interior temperature, and r is half the total vertical 
optical depth. The inner parts of black hole accretion disks have opacities 
that are generally dominated by electron scattering, so r ~ ktS/ 2, where 
E is the surface density of the disk, which is constant on these time scales 
very much less than the radial flow time scale. The Thomson opacity kt is 
also constant, being independent of temperature provided there is already 
sufficient ionization (and there is), so the optical depth is independent of 
temperature and the cooling rate per unit area is F~ oc T 4 . Note that I 
have made the usual assumption here that F~ is a real heat flux and does 
not involve other forms of energy. 

The dissipation rate per unit area is 


F + 


T dn 


(2.5) 
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Vertical hydrostatic equ ilibr ium implies that the disk half thickness H 
2P/(HjlE), so that eq. (2.5) becomes 


F+ rv 2RP T R<I> 


D|E d R' 


( 2 . 6 ) 


In the radiation pressure dominated inner region, P ~ aT 4 /3, so that equa¬ 
tions (2.6) and 0 ) or (^3|) imply that F + oc T 8 ! Hence a perturbative 
increase in temperature increases both the local cooling and heating rates, 
but the heating rate increases much faster, leading to a thermal runaway. 

Such thermal (and “viscous”) instabilities arising from opacity varia¬ 
tions in hydrogen ionization zones rather than radiation pressure have been 
applied with considerable success to understanding the outburst behavior 
of dwarf novae (e.g. ^9|) and soft X-ray transients (e.g. |^]). The reality of 
these instabilities in radiation-pressure dominated zones has never been es¬ 
tablished, however. It has been argued j9^] that the stress should be propor¬ 
tional to the gas pressure only, in which case r ^ oc pT oc E 2 f2|-T/(4P) oc 
T -3 in the radiation pressure dominated inner zone. This implies that F + 
depends only linearly on temperature, producing a thermally (and in fact 
“viscously” as well) stable flow. We will examine this argument in section 
5 below. In addition, we will look at how radiation magnetohydrodynamic 
instabilities might affect thermal instabilities in accretion flows. 

Observationally motivated models of accretion flows generally start with 
the assumption that the flow is stationary and axisymmetric about the ro¬ 
tation axis, even though whatever is responsible for the anomalous angular 
momentum transport must almost certainly involve time-dependent, non- 
axisymmetric fluctuations. (Time-dependent models are often constructed 
to model the response of the disk to thermal and “viscous” instabilities.) 
Another simplifying assumption is that vertically-integrated flow equations 
provide a reasonably accurate description of the behavior of conserved quan¬ 
tities, even if the flow is not very geometrically thin. The resulting steady- 
state conservation laws are those for mass, 


M = AnRHpv , 


(2.7) 


radial momentum, 




( 2 . 8 ) 


angular momentum, 




(2.9) 
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and energy, 


- 4nRF adv = M 


AU 

dR 


P- 


AR 


= 4t tR 2 Ht r ^ + 2nR(2F~) 
art 


( 2 . 10 ) 

Here M is the constant accretion rate through the flow, R is the distance 
from the rotation axis, v(R) is the inward radial flow speed, £{R) is the 
angular momentum per unit mass, U(R) is the internal energy per unit 
mass, and F~(R) is the energy flux leaving the flow on each of the two 
vertical surfaces. 

Note that the anomalous stress T R( j, is assumed to enter the angular 
momentum and energy equations in exactly the same way as an ordinary 
fluid viscous stress: torques are exerted in a dissipative fashion. This need 
not be the case. External torques exerted on the disk by a global, ordered 
magnetic field in a magnetohydrodynamical (MHD) wind can in principle 
extract energy and angular momentum in precisely the right ratio necessary 
to allow material to accrete with no dissipation |n]]. The disk can be ice 
cold and still accreting! Because of possibilities like this, I have been careful 
to call F~ an energy flux and not a heat flux. Even if internal turbulence 
is responsible for angular momentum transport, it is still not clear that all 
the accretion power should be converted into heat. Instead, some of that 
power may be converted into bulk kinetic and magnetic energy, and F~ 
must include these non-radiative contributions. 

Equations (2.S) and ( 2.10| ) diffe r from those of standard geometrically 
thin accretion disk theory (e.g. [ 102]) by the inclusion of advective and radial 
pressure support terms which have received huge theoretical attention in 
recent years in the guise of advection dominated accretion flows (ADAFs, 
Q, J|(J). In the absence of these effects, one gets the standard equations of 
thin disk theory, i.e. that the angular velocity is Keplerian, 


Si(R) = fi K (R), 


( 2 . 11 ) 


and local internal dissipation is balanced entirely by local vertical cooling, 

dflx 


F = -RF[T R(t> ■ 


AR 


( 2 . 12 ) 


The angular momentum equation (T9) can be integrated from an inner 
radius of the flow R m (usually taken to be the innermost stable circular 
orbit, or ISCO, radius for black hole accretion disks) out to an arbitrary 
radius R to give 


M{1~ 4i) = 4?r R 2 Ht R 4 , - ATTR? n H in T R(/> u 


(2.13) 


Combining this with the local energy balance condition ( [2. 12 ) for geometri¬ 
cally thin disks allows us to write down an expression for the emergent flux 
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which is independent of the anomalous stress throughout most of the flow, 


F~ = 


3GMM 

8nR 3 


1 - 


Rm 

R 


1/2 


47T R ]u H' ln T S'^,in 

M{GMR)V 2 


(2.14) 


This is the most beautiful result of standard accretion disk theory: that 
steady state accretion disks have an emergent flux that is independent of 
the details of the anomalous stress. Of course, accretion flows around black 
holes are observed to vary on all sorts of time scales, so steady state may 
not be a good assumption. Moreover, there is still a dependence on the 
stress at the inner edge of the flow. Until recently, this stress was 

almost always assumed to vanish near the ISCO, due to the fact that the 
inward transonic flow inside this point was presumed to rapidly decouple 
from the disk. This assumption has been challenged recently , and the 

physics of angular momentum transport near the ISCO has since received 
considerable theoretical scrutiny. We will examine this further in section 4 
below, and for now, we will keep the inner stress in our equations. 


Equation (2.14) can be integrated to give the total luminosity of the 


disk. 


L = 


r°° GMM 

/ 2nR(2F~)dR = —— + (4^ rl H in T R0m )Q m . (2.15) 

J Rin 


This has a simple physical interpretation. The first term is the binding 
energy per unit mass of material at the inner radius, times the accretion 
rate. The second is the rate at which work is being done on the disk at 
the inner radius. By defining the accretion efficiency in the usual way, 
rj = I/(Mc 2 ), with i]ss being the efficiency under the standard (Shakura- 
Sunyaev) assumption of a no-torque inner boundary condition and Arj = 
rj — riss being the additional efficiency due to the inner torque, equation 
(2.14) can be written in a physically appealing way (e.g. [7|), 


F~ = 


3 GMM 
87 tR 3 


Rr 


1/2 


'-'if +4 nf i 


i?i, 


Ru 


1 / 2 ' 


(2.16) 


Here r g = GM/c 2 is the gravitational radius. (The full relativistic version of 


eq. [p,16| may be found in [IT].) Note that the effects of the inner boundary 
on the flux decay with radius as i?~ 7 / 2 , only slightly faster than the overall 
asymptotic decay of the flux at large radius, oc R~ 3 . 

So far we have been neglecting the advective terms in our conservation 
laws ( |2.7| )-( ^T0 ). If we include them, equation ( 2.14 ) becomes 


F~= — 


Ml dD 
47rl? di? 


1 - -=■ 


fin 4?r R ]n Il Ui T Ro ^ 


Ml 


Eadv- 


(2.17) 
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Advection, which becomes increasingly important as one moves inward in 
the flow, has two effects: the rotation profile is no longer Keplerian, and the 
orbital energy dissipated at any radius is partly advected (lost!) r adiall y 
inward. The latter effect is evident from the last term in equation (2.17). 
The emerging energy flux F~ is less than it would be if advection was not 
included. 

For luminous accretion disks, advection starts to play a big role when 
the overall luminosity starts to get above some fraction (~ 0.3 in standard 
models) of the Eddington luminosity LEdd- The flow is then geometri¬ 
cally thicker, becoming a “slim” or “thick” accretion flow. Such flows are 
thermally and “viscously” stable [il|.^], because F a d v has a temperature de¬ 
pendence strong enough to beat the temperature dependence of F + . The 
radial and vertical structure of these flows is no longer decoupled, and the 
emergent flux F~ now depends on the very uncertain assumptions that are 
usually made to solve for the vertical structure. For a discussion of the 
fully relativistic equations of advective accretion flows, see e.g. ||, ||, [ p2| , 
and j|(|. Fully relativistic models of optically thick, “slim” accretion disks 
have been constructed in [^?) . Models of hyper-Eddington accretion flows 
onto stellar mass black holes relevant for gamma-ray burst models have also 
been constructed [Q. By far and away most of the energy dissipated in 
such flows is advected into the black hole until the accretion rate exceeds 
~ 10~ 2 M 0 s , where neutrino losses can cool the gas on the flow time 
scale. 

Before leaving this section, I should mention the considerable atten¬ 
tion that has been devoted recently to optically thin, radiatively inefficient 
ADAF solutions around black holes. Because we are focusing on luminous 
accretion flows here, I will have little to say about these solutions, but it is 
important to recognize that they may nevertheless arise over certain ranges 
of radii even in otherwise luminous flows. In particular, it is noteworthy 
that such flow solutions exist that are thermally and “viscously” stable over 
much the same range of accretion rates that geometrically thin, optically 
thick disk solutions exist (e.g. Q), and it is not clear why nature should 
choose one flow solution over the other. The presence or absence of an inner 
ADAF flow might explain the different spectral states observed in black hole 
X-ray binaries J39|. One should bear in mind, however, that all the models 
discussed in this section are based on the a stress prescription, and it is fair 
to say that the microphysics of optically thin ADAFs is even more poorly 
understood than that of optically thick flows (e.g. [^5), |k]]). 


3 Spectral Formation 

As we saw in the last section, if we neglect the effects of advection, then 
the energy flux emerging from each face of a stationary accretion disk can 
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be written 

, , 3 GMM , i s , s 

F ~ (R) = (■ R/r g ), (3.1) 

where T a i^i(x) is a function which in full general relativistic treatments de¬ 
pends on the spin parameter a/M of the hole but still approaches unity at 
large radii. In standard models with a no-torque inner boundary condition, 
X a /M a l so vanishes at the ISCO, but as we have seen, this need not be the 
case. Once again, let me remind the reader that although F~ is usually 
assumed to be the emerging radiation flux from the disk surface, this need 
not be the case. It is actually the total energy flux, and may include contri¬ 
butions from magnetic and bulk kinetic energy. We will return to this point 
shortly, but for now, let us make the standard assumption that F~ really 
does represent an emerging radiation flux. 

What, then, is the emerging spectral energy distribution? An observer 
here on earth, viewing the disk at a distance d away at an inclination angle 
i to the rotation axis, will measure a spectral energy flux 

COS 1 f ^ oul 

F u = —p-J 2nRI v (R,i)dR, ( 3 . 2 ) 


where /„ is the angle (limb darkened or brightened!) and radius-dependent 
local emergent specific intensity, and we have introduced a finite outer ra¬ 
dius I?out to the disk. (Once again, we are using Newtonian physics here, 
and neglecting e.g. Doppler shifts, gravitational redshifts, the gravitational 
bending of light rays, and the relativistic proper emitting area, but these 
are all straightforward to account for using relativistic “transfer functions”, 
e.g. | 

The simplest (and therefore most often used!) way of estimating the 
shape of the spectrum is to just naively assume that each annulus radiates 
like an isotropic local blackbody, i.e. I v (R,i) = B v [T e (R)\, where B v is 
the Planck function with temperature given by the effective temperature 
determined from equation ( |.3.l| ), T e (R) = [F (R)/^] 1 / 4 . If, say, the radial 
dependence of the emerging radiation flux can be approximated as a power 
law, F~(R) cx R-P, then a simple change of integration variable to x = 
hv/(kT) implies that the shape of the observed spectrum is given by 


») 


F v oc v 3 - 8 ' 13 


nhv /(fcTeout) 


' hv/(kT e in) 


x S/0 l^ x 
e x - 1 


(3-3) 


As one would expect, at low frequencies satisfying hv <!C kT eout we simply 
end up with F v cx v 2 (the Rayleigh-Jeans portion of the spectrum emitted by 
the outermost annulus), while at high frequencies satisfying hv kT e i n we 
have F v tx v 3 exp[— hv/(kT e i n )\ (the Wien portion of the spectrum emitted 
by the innermost annulus). In between, for kT eout <C hv <C kT e j n , the 
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Fig. 1. Illumination geometry for a source on the symmetry axis at height H * 
above the equatorial plane. The angle a between the disk normal and the incoming 
light ray is simply 7 t/2 — (/3 + 9), where 9 is the angle between the tangent to 
the local disk surface and the equatorial plane, and /3 is the angle between the 
illuminating light ray and the equatorial plane. For large R, /3 and 9 are both 
small angles, with /3 ~ (//* — H ) /R and 9 ~ AH/AR. 


observed spectrum is a superposition of blackbodies from many different 
radii, producing a power-law spectrum F v oc . Now, at large radii, 

equation ( [j.l| ) implies that P = 3, giving the famous result that F v oc v 1 / 3 
for a multi-temperature blackbody accretion disk. 

There are many ways in which this standard result can be modified. One 
possibility is that energy liberated at small radii in the flow is reprocessed 
by the disk at larger radii. By way of illustration, consider a source of 
luminosity L* located on the rotation axis at height iJ* above the equatorial 
plane. Figure [j] illustrates the geometry of the problem. 

Neglecting general relativistic effects, in particular light bending, the 
irradiating flux on the outer disk at radius R will be given by 


F irr = L H R2 a) COS a, (3.4) 

where a is the albedo of the outer disk, averaged over the illuminating 
spectrum, and a is the angle between the disk normal and the incoming 
light ray. Taking the illuminating radius R to be much larger than the 
local disk height H , the source height iJ*, and also dH/dlnR, then a little 
trigonometry gives 


_ L*(l-o) (H\ (dlnH i?*\ 
4 t tR 2 \R ) \d\n R + H ) ' 


(3.5) 


For a source located near the equatorial plane, e.g. the surface of the 
inner disk itself, we have ii* <C H(R). Equation (3.5) then implies that 
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significant reprocessing will only occur if the outer disk geometry flares up, 
with H/R an increasing function of R , or is warped |m|]. In this case, pro¬ 
portionately more power is transferred to larger radii, and this decreases the 
value of (3 and reddens the spectrum. Standard models of disks work in just 
this way, with H <x R 21 / 20 and R 9 / 8 for the so-called “middle” and “outer” 
regions, respectively. This gives an irradiating flux approximately propor¬ 
tional to R~ 2 , which must therefore dominate the more rapidly decaying 
(oc R ~ 3 ) local dissipation, and giving F v oc v~ x (e.g. ]4^1). 


For a sour ce located far above the equatorial plane, with H * H(R ), 
equation (3J5) gives an irradiating flux oc R~ 3 , producing the same F v oc 
v 1 / 3 spectrum as a standard disk but with enhanced power. One way in 
which this might be done is by having a geometrically thick, hot corona or 
ADAF in the inner parts of the flow. General relativistic light bending can 
also produce a similar effect: an observer located on the outer disk can see 
the opposite side of the disk at an effectively high if* [Q. This effect alone 
can produce significant additions to the locally dissipated power in the case 
of rapidly rotating Kerr black holes [B5|. 


Even without reprocessing, the emergent spectrum can be affected by 
alterations in the radial distribution of the local radiative flux F~. For 
example, if an inner torque across the ISCO produces a large increase in the 
accretion efficiency, equation ( |2.16 ) implies that a substantial range of inner 
radii could have F~ oc R~‘^ 2 , giving a spectral energy distribution from 
equation ( |i..i| ) of the form F v oc 1 / 5 / 7 , bluer than i/ 1 / 3 ! On the other hand, 
gravitational light bending from the much brighter inner disk enhances the 
amount of reprocessed flux at larger radii, giving F] rr oc R~ 3 and enforcing 
a more usual F v oc v 1 / 3 spectrum at longer wavelengths |7j]. 

Advection also alters the radial distribution of F ~, causing it to rise less 
steeply toward smaller radii because more of the accretion power is advected 
into the black hole rather than being radiated. This in turn implies that 
F~ falls less steeply when moving out in radius, implying a redder spectrum 
than F v oc v 1//3 [ 1151. In fact, models of slim disks can give flux distributions 
as flat as F~ oc R~ 2 , and equation ( fUi| ) then gives F v oc u~ l [121,122, Bl[. 
The global effect of advection on the spectrum is that the luminosity is 
less than you would expect for the given accretion rate, because advection 
necessarily reduces the radiative efficiency. 


Most models of the vertical structure of accretion disks around black 
holes imply that local blackbody emission is likely to be a very poor ap¬ 
proximation. The most serious problem with it is the neglect of electron 
scattering, which is often far greater than true thermal absorption opacity. 
Again, a crude approach to handli ng th is has been used from the very be¬ 
ginning of accretion disk modeling [102], and that is to use a so-called local 
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modified blackbody spectrum at each radius, 


where 


r 2 B„{T)el /2 
" 1 + ej 2 ’ 

(3.6) 

th 

^ A-v 

e v — 1 til 

n e a T + Xl 

(3.7) 


is the photon destruction probability, with 
coefficient, n e the electron number density, 
section. 


being the thermal absorption 
and <jt the Thomson cross- 


The physics behind equation (3.6) lies in the idea that, once created by 
thermal emission processes, photons random walk out of the atmosphere. If 
absorption opacity is smaller than electron scattering opacity, this random 
walk does not destroy photons generated within an effective optical depth 
"Jeff = (TabsTT) 1 / 2 of the surface, where r a b s is the true absorption optical 
depth and tt is the Thomson depth. Assuming LTE and integrating the 
emission over this vertical layer immediately gives I v ~ B u ef , in rough 
agreement with equation (3.6). Provided the accretion disk is effectively 
thick (r e ff > 1 at the frequency of interest), then the spectrum may be 
viewed as coming from a depth in the disk corresponding to r e s = 1. As 
shown in Figure the innermost regions of some black hole accretion disk 
models are not even effectively thick, however, and the observed spectrum is 
formed throughout the entire geometrical thickness of the disk. In such cases 
the spectrum will be directly sensitive to the assumptions that went into 
creating the entire vertical structure. Even when r e g > 1, the spectrum will 
depend on the ambient densities and temperatures at the r e fj = 1 effective 
photosphere, and these will in turn depend on the vertical structure of the 
disk model. 

Note that the modified blackbody spectrum in equation (3.6) has I v < 
B V {T). This makes sense, as deviations from blackbody imply that the 
disk is a thermodynamically less efficient radiator. Because it must still 
radiate the same total flux, the ambient gas temperature must be hotter 
than the effective temperature, and the resulting spectrum therefore extends 
to higher photon energies. When integrated over the disk, the spectral 
energy distribution must therefore flatten below the canonical F v oc v 1 / 3 as 
photons are redistributed to higher energies. These effects are illustrated 
in Figure |^, which compares stellar atmosphere calculations of relativistic 
disk spectra with the local blackbody assumption. 

In addition to the effects of Thomson scattering on the spectrum, Comp¬ 
ton scattering can also be important in hot disks. The net Compton heating 
rate of the plasma per unit volume is 

hv 4 kT \ / 47tJ n 


n e <7TC 


m e c“ 


m e c 




(3.8) 
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Fig. 2. The midplane effective optical depth as a function of radius in relativistic 
accretion disk models around a ten solar mass black hole with L = 0.3Z/Edd and 
TR'j, = aPtot, where Ptot is the total (gas plus radiation) pressure. Both Kerr 
( a/M = 0.998) and Schwarzschild models are shown, and a no-torque boundary 
condition is assumed at the ISCO, which is why the curves rise upward at small 
radii. The effective optical depth is actually a function of frequency, and what is 
depicted here is an average using the Rosseland mean free-free absorption opac¬ 
ity. Solid and dashed curves assume a = 0.1 and 0.01, respectively. The lower 
anomalous stress of the latter implies higher surface density and higher effective 
optical depth in both the Kerr and Schwarzschild cases. 


where J is the frequency-integrated mean intensity of the radiation field, 


nOO -j^ rOC /» 

J= dvJv = — / dv ® 

Jo 4 tt Jo J 


dttlv, 


and v is a specially defined average frequency, 

/•OO 

to 


1 

V = — 

J. 


^d + ^). 


(3.9) 


(3.10) 


(The second term in parentheses in equation | 3.1C | represents the effects of 
stimulated scattering.) The Compton heating rate can be compared with 
the net heating rate from thermal absorption and emission. Assuming LTE, 
this is 

t) 


fOO 


th 


c / d,vx'J 
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— B„). 


(3.11) 
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Fig. 3. Spectral energy distributions of relativistic accretion disk models around 
a 1.25 x 10 8 Mq Kerr ( a/M = 0.998) hole, accreting at 0.25 Mg yr _1 , viewed at 
60° from the rotation axis J(i^]. The dashed curve, which extends to the highest 
energies, assumes that electron scattering is coherent (Thomson scattering). The 
solid curve allows for the effects of Compton scattering. Compton downscattering 
pulls the high energy tail of the spectrum down, while Compton upscattering 
increases the spectral luminosity at lower energies. The dot-dashed curve is the 
spectrum that results from assuming blackbody emission at the local effective 
temperature at every radius, and is clearly a very poor approximation. 


Defining the Thomson opacity as 


_ n e o t 

Avt = - 

p 

and the Planck mean opacity as 


ftp = 



(3.12) 


(3.13) 


we see on comparing equations (3.8) and ( |3.11 ) that Compton scattering 
will dominate the thermal coupling between gas and radiation if 


ftT 


4 kT 

m e c 2 


> 


Kp. 


If, in addition, the effective Compton y-parameter, 

4 kT t 2 \ 

Veff = -jnrax (r T , t t ) , 

m e c* 


(3.14) 


(3.15) 
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Fig. 4. Spectral energy distributions of accretion disk models around Kerr ( a/M = 
0.998) black holes at fixed Eddington ratio viewed nearly face on, for various 
black hole masses: 10 s (top), 10' (middle), and 10 6 Mq (bottom) Jb^]. Note the 
increased relativistic smearing of the Balmer and Lyman edges of hydrogen for 
higher black hole mass. This is because the temperature of the disk is lower at 
higher mass, so these features are formed closer to the black hole. 


exceeds unity, at the particular frequency of interest, then Compton scat¬ 
tering will modify the spectrum, either by downscattering if hv > 4fcT, or 
upscattering if hv < 4 kT. (In equation 1 3.15 ], tt is evaluated at the T e ff = 1 
surface.) Figure || illustrates these effects on the spectrum. 


Numerous authors have attempted detailed calculations of the spec¬ 
tral energy distributions emerging from accretion disks around black holes 
using stellar atmosphere modeling at every radius of the disk (see e.g. 
|pd^,|7l|,p7|,|l0^, pT2| , |l07| ,[37|,|63|,[l0l|,[64|,|6l|,|6^| for models in the AGN con¬ 
text). Many of these models fully incorporate general relativistic effects on 
the disk structure, and also include relativistic Doppler shifts, gravitational 
redshifts, and gravitational bending of light rays. These latter effects play an 
important role in smearing out atomic absorption and/or emission features 
in the spectrum, as illustrated in Figure f|. Many models also incorporate 
sophisticated treatments of the radiative transfer, Compton scattering, and 
non-LTE effects in the atomic level populations of numerous elements and 
ions. 


As sophisticated as these models are, however, they still rest on variants 
of the Shakura-Sunyaev anomalous stress prescription, and are subject to all 
the uncertainties we discussed above in section 2. Moreover, just as a stellar 
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model must specify a distribution of nuclear energy generation, a model of 
the vertical structure of an accretion disk must specify the distribution of 
turbulent heating. The presence of turbulence in the disk also begs the 
question as to whether heat is transported radiatively, or whether bulk 
transport in the turbulence itself also plays a role. As we discuss below in 
section 5, standard assumptions often used to model the inner, radiation 
pressure dominated regions lead to equilibria that are convectively unstable 
|p8|| , although it is far from clear how that convection would manifest itself 
in a flow that already requires turbulent angular momentum transport. 


Perhaps the most serious question plaguing current Shakura-Sunyaev 
based spectral models is the partitioning of the dissipation within the op¬ 
tically thick and thin regions of the disk. There is widespread empirical 
evidence for the existence of Comptonizing hot plasma in both black hole 
X-ray binaries and AGN, and one possibility is that this plasma exists as 
a magnetized corona above the disk photosphere. After all, the turbulent 
convection zone of our own sun generates a hot corona, so why shouldn’t 
a turbulent accretion disk have one as well? (This is one version of in¬ 
cluding non-radiative energy in F ~.) Models have even been considered in 
which most of the turbulent dissipation is assumed to occur in the corona, 
and not in the disk interior j49|! If a substantial fraction of the accretion 
power is in fact dissipated outside the disk, then there are numerous con¬ 
sequences 1114 : the disk becomes geometrically thinner and denser, more 


effectively optically thick, and more gas pressure dominated. 


Accretion disk physics is uncertain enough, but the physics of the hot 
plasma that must exist in real sources is even more insecure. Whether 
or not the hard X-rays observed in black hole X-ray binaries or AGN are 
generated in a disk corona or elsewhere in the flow (e.g. an inner ADAF that 
we discussed at the end of section 2), some of these X-rays will illuminate 
the disk from the outside, and thereby be reprocessed. 


One of the most interesting X-ray reprocessing features is fluorescent 
iron Ka line emission. The great excitement surrounding this line is that 
in at least three sources (MCG-6-30-15 [116, 123| , |jj, Mrk 766 [f78| , and 
NGC 3516 ||84| , |l2C|| ), it is clearly observed to be relativistically broadened. 
The resulting line profile is a convolution of the spatial emissivity profile of 
the line with relativistic Doppler shifts, gravitational redshifts, and grav¬ 
itational light bending |70|. In other words, the line profile gives us the 
unprecedented opportunity to map out the equatorial test-particle orbit 
structure in a black hole spacetime, assuming the emitting material (the 
disk) is geometrically thin. 

The Ka line is produced because hard X-ray photons with energies above 
the K-edge of iron can knock inner shell electrons out of the atom. The 
resulting ion is born in an excited state and can return to the ground state by 
having an L-shell electron drop down to the vacancy in the K-shell, emitting 
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either a Ka line photon or one or more Auger electrons. The strength of 
the resulting emission line is sensitive to the ionization state of iron in the 
outer layers of the disk. At low ionization levels, where both the K and L 
shells of most iron ions have no vacancies, the material is transparent to Ka 
photons which then escape more or less freely from the disk once they have 
been created by X-ray fluorescence. Once iron is sufficiently ionized that 
vacancies in the L-shell start to appear (FeXVIII to FeXXIV), then a Ka 
photon can resonantly scatter off the iron atoms, enhancing the probability 
of its destruction through the Auger process. The emerging line strength is 
therefore reduced. This continues until iron is sufficiently ionized that only 
FeXXV and higher is present. The Ka line photons still resonantly scatter 
with such ions, but the Auger effect cannot happen because there are no 
L-shell electrons to be ejected. At very high ionization when iron is fully 
stripped and the recombination rate is small because of high temperatures, 
no photoionization can occur and the emission line disappears. For an 
excellent review of the physics of the iron Ka line, see p(i| |. 

The energy powering the iron Ka line emission comes from photoelectric 
absorption of X-rays above the iron K-edge. Such photoelectric absorption 
by iron and other heavy elements present in the disk is in fact an efficient 
means of converting illuminating X-ray power below ~ 10 keV into heat 
in the outermost layers of the disk Figure || shows the albedo 

of an illuminated slab of material as a function of incident photon energy 
jnj. Below ~ 10 keV, most of the incident X-rays are absorbed due to 
photoionization of heavy elements, and the level of that absorption depends 
on the ionization state of the disk. Above ~ 10 keV, the albedo quickly rises 
to near unity as the incident X-rays are simply reflected back by Thomson 
scattering. This is widely believed to explain the characteristic upturn of 
the X-ray spectrum that is observed at these energies in both AGN and 
black hole X-ray binaries | 
of relatively cold material in these sources. (Note that especially in the case 
of AGN, however, this cold material need not be the disk, but instead could 
be gas much further out from the source.) At higher energies the albedo is 
again reduced to well below unity, independent of the ionization state of the 
slab, because X-rays are Compton downscattered and therefore lose energy. 

Just as the thermal emission spectra from accretion disks can depend 
sensitively on the details of the vertical disk structure, X-ray reprocess¬ 
ing and reflection features also depend on the assumptions one makes in 
describing the physics of the illuminated slab. Models that incorporate 
hydrostatic equilibrium produce surface temperature profiles containing a 
number of sharp discontinuities due to transitions between different regions 
of thermally stable material |87],|22|], and the resulting reflection/fluorescent 
spectra are altered by this fact. 

In addition to the very uncertain geometry of the illumination itself, one 


5,124], and is strong evidence for the presence 




20 


The title will be set by the publisher. 



Fig. 5. Energy-dependent albedo of an X-ray illuminated, cold slab with solar 
abundances, viewed at twenty degrees from the slab normal. The solid curve 
assumes the slab is neutral, apart from hydrogen and helium which are fully 
ionized. The dashed curve is for an ionized slab, with ionization parameter £ = 
100. In both cases the illumination was assumed to be by a power-law spectrum 
with photon index -2 and exponential cutoff energy of 100 keV. (This figure was 
kindly generated by Shane Davis, courtesy of the XSPEC routines pexrav and 

PEXRIV @0 


wonders whether even this treatment of the physics is enough. The bottom 
line of this section is that many aspects of our models of observed accreting 
black hole sources are fundamentally limited by uncertainties in the basic 
physics of the accretion flow itself. For the rest of these lectures we will 
switch gears and address this physics itself, explaining how to go beyond 
Shakura-Sunyaev type models. 

4 The Physics of Angular Momentum Transport 

A system with fixed total angular momentum in thermal equilibrium must 
be in a state of rigid body rotation This fact lies at the very heart of 
the way accretion disks work, because it shows that differential rotation is 
a source of thermodynamic free energy in the flow. As discussed long ago 
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by Lynden-Bell & Pringle j7f|, one can always reduce the total rotational 
kinetic energy of an accretion disk by fixing its mass distribution and redis¬ 
tributing the angular momentum so that it is uniformly rotating. However, 
such a state is generally incompatible with mechanical equilibrium in the 
gravitational field of the central object, which requires that the angular ve¬ 
locity decrease outward. In this case, pairs of orbiting fluid elements moving 
with different angular velocities can lower their mechanical energy by shift¬ 
ing angular momentum from the element with higher angular velocity to the 
element with lower angular velocity |f75f . This transports angular momen¬ 
tum outward through the disk, causing material to lose rotational support 
and flow in to regions of greater binding energy. This general behavior 
is simply the action of the second law of thermodynamics, and is entirely 
analogous to the one-way flow of heat from regions of high temperature to 
low temperature. Ordinary microscopic viscosity acts to transport angular 
momentum in exactly this fashion, but this is not sufficient for astrophysical 
accretion disks. 

Some form of non-microscopic, “anomalous” angular momentum trans¬ 
port mechanism or mechanisms must therefore be at work, and it is impor¬ 
tant to recognize that it need not be the same in all accretion flows in the 
universe. Inherent in the Shakura-Sunyaev “guess” as to the form of the 
anomalous stress T r( j, is the existence of some form of turbulence in the flow: 
ideally a form of turbulence that is generated by the differential rotation 
itself. Turbulence is often, though not always, generated when a laminar 
flow is dynamically unstable to perturbations. These perturbations may be 
infinitesimal, so that the slightest touch of a feather will produce instability. 
The initial growth of the instability may then be understood from linearized 
equations about the equilibrium flow, which is then said to be linearly unsta¬ 
ble. Linear equations are much easier to solve than nonlinear equations, and 
much of the theoretical work has therefore been devoted to linear instabili¬ 
ties. Understanding the response of the flow to nonlinear, finite amplitude 
perturbations is a lot harder, but there do exist laminar flows that are non- 
linearly unstable even while being linearly stable. Of course to theoretically 
understand how a flow really changes in response to the development of even 
a linear instability requires a treatment of the nonlinear equations, which 
in practice is usually accomplished through numerical simulation. 

The simplest way to examine the linear stability of a fluid is to infinitesi¬ 
mally perturb a fluid element locally and examine whether the forces exerted 
on the element act to return it to its original position. Stability criteria de¬ 
rived in this way are necessary, but not sufficient because it is also possible 
for the global structure of the flow to produce instabilities. Unfortunately, 
analyzing the global stability of an accretion flow is technically much more 
difficult, even within linear theory. One must solve a global eigenvalue 
problem for the linear wave modes of the system, fully accounting for the 
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boundary conditions on the flow. In order to derive a stability criterion, 
one must also ensure that the modes examined form a complete set so that 
any initial perturbation data can be represented in terms of them. 

Global stability analyses generally recover the results of local stability 
analyses in the limit of very short wavelength perturbations (often called 
the WKB limit in the literature). If global unstable wave modes are found, 
it may be that changes in the assumed boundary conditions can stabilize 
the modes. This can never be true of local instabilities, however, which care 
only about local conditions in the flow. In this sense local instabilities are 
much more robust and generic than truly global instabilities. 

In my discussion of accretion disk instabilities in this section, I will 
always neglect the self-gravity of the accreting plasma itself. Accretion rates 
in both AGN and black hole X-ray binaries are thought to be sufficiently 
small that the tidal field of the black hole completely overwhelms the self¬ 
gravity of the disk, at least near the hole where most of the luminosity is 
generated. Self-gravity is however more important further out in the disk 
in AGN, and may also be important in hyper-Eddington flows invoked in 
gamma-ray burst models. 

Within ideal hydrodynamics, local linear stability of an axisy mmet ric 
rotating flow is guaranteed if the Hpiland criterion is satisfied (e.g. EH): 



N 2 + k 2 >0 

(4.1) 

and 

{VP x R) • (W 2 x VS) > 0, 

(4.2) 

where 

2 1 dl 2 1 ~ - 

K ~ wlm ~ rz r ' 

(4.3) 

is the square of the epicyclic frequency and 



N 2 = -VP ■ (-Vp - —VP ] = —VP ■ vs 

P \P l p ) pep 

(4.4) 


is the square of the Brunt-Vaisala frequency. Other quantities in these 
equations are the pressure P, the density p, the entropy per unit mass S, 
the heat capacity per unit mass at constant pressure cp, the ratio of heat 
capacities 7 , the angular momentum per unit mass l, the distance from the 
rotation axis P, and the corresponding cylindrical polar coordinate radial 
unit vector R. (A relativistic generalization of the Hpiland criterion can be 
found in 1001 .) 

The Hpiland criterion is easy to understand in two limits. For non¬ 
rotating equilibria (e.g. a non-rotating star), the criterion reduces to the 
Schwarzschild criterion that the entropy must not increase inward for sta¬ 
bility against convection. Or, perhaps in more familiar form, the tempera¬ 
ture must not increase inward faster than the adiabatic (constant entropy!) 
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temperature gradient. Provided this is true, local fluid elements will simply 
oscillate at the local Brunt-Vaisala frequency under stable buoyancy forces, 
and these restoring forces are responsible for supporting gravity waves (g- 
rnodes) in stars. To isolate the effects of rotation, consider an equilibrium 
that has constant entropy everywhere. Then the Hpiland criterion reduces 
to the Rayleigh criterion: the specific angular momentum must not decrease 
outward for stability. Physically, if one perturbs a fluid element radially 
outward, it conserves its own specific angular momentum. If the ambient 
specific angular momentum decreases outward, then the fluid element will 
be rotating too fast to stay in its new position, and centrifugal forces will 
fling it further outward. Stability in this case implies that perturbed fluid 
elements will oscillate back and forth at the local epicyclic frequency, and 
organized epicyclic oscillations form the basis for Rossby waves (r-modes) 
in rotating stars. It is no accident that g-mode and r-mode physics appear 
together in the Hpiland criterion, as these modes are fundamentally inter¬ 
twined when both angular momentum and entropy gradients are present. 

The Hpiland criterion is a huge disappointment for understanding why 
turbulence might exist in accretion disks. The piece of it that is related 
to the rotation is connected to specific angular momentum gradients, and 
these are strongly stable for rotation profiles that do not differ too much 
from Keplerian. Of course, the Hpiland criterion is only a local stability 
criterion. There do exist global, linear hydrodynamic instabilities in ac¬ 
cretion disks. The one I am most familiar with is the Papaloizou-Pringle 
instability |)lj but, like all global instabilities, its existence is sensitive to 
the assumed boundary conditions p9| . Its nonlinear development also ap¬ 
pears to produce large scale spiral waves rather than local turbulence [ ^o[ , 
although such waves are perfectly capable of producing “anomalous” angu¬ 
lar momentum transport. The Hpiland criterion also only addresses linear 
stability. Hydrodynamic, planar shear flows can be violently unstable to 
finite amplitude perturbations even if they are linearly stable, and this sug¬ 
gests that something similar might happen in accretion disks. However, 
extensive searches for such instabilities with numerical simulation have con¬ 
sistently failed to find them in linearly stable flows, except very close to the 
limit of marginal stability |^o|,|5^]. It is noteworthy that this is true even 
though the same simulation hardware easily and correctly finds the insta¬ 
bilities in planar shear flows J2(], ^6). As a point of principle, it is possible 
that the simulations are somehow not accessing an unknown destabilizing 
mechanism (see |7j| for a recent argument in this direction). However, one 
should keep in mind that there is an important physical difference between 
linearly stable planar shear flows and differentially rotating flows: the latter 
have a strong local restoring force arising from specific angular momentum 
gradients |f 2 (j|] . 

There is a key piece of physics that is central to the Hpiland criterion: 
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perturbed fluid elements conserve their own entropy and specific angular 
momentum. This is not true when we include the effects of magnetic fields, 
as they fundamentally alter the local dynamics of perturbed fluid elements. 
To order of magnitude, the characteristic frequency associated with mag¬ 
netic fields, within the MHD approximation anyway, is simply kv a, where k 
is the perturbation wavenumber and v\ = B/(Airp) 1 / 2 is the Alfven speed. 
Magnetic fields will alter the character of hydrodynamic r- and g-modes 
whenever kv a is of order the relevant frequency (|k| or \N\) of these modes. 
(This is true even when the hydrodynamic modes are unstable, e.g. mag¬ 
netic fields can strongly affect hydrodynamic convective instabilities which 
have imaginary N.) At first sight one would think that a weak magnetic 
field (low Alfven speed) would therefore have no effect. However, no matter 
how weak the magnetic field is, one can always make it important by going 
to sufficiently high wavenumbers (short length scales), e.g. k ~ |k|/va in 
the case of r-modes. Indeed, for k |k|/ua, one expects rotation to be 
irrelevant as the modes become magnetically dominated, no matter how 
weak is the magnetic field. This all assumes ideal MHD, in particular that 
magnetic field lines are frozen into the fluid. Resistive effects, which allow 
the field lines to “slip” with respect to the fluid, become more important on 
small scales. If the field is so weak that it only becomes important below 
resistive length scales, then we are back to hydrodynamics on all scales. 

The basic physics as to how weak magnetic fields destabilize accretion 
disks is most transparent using a mechanical analogy developed by Balbus & 
Hawley |l7| . Consider two identical test particles of mass m moving together 
(right on top of each other) in a circular orbit in an external gravitational 
potential <f>(.R). The angular velocity H of the particles’ motion will of 
course be given by 


1 <9<f> 
RdR’ 


(4.5) 


Now suppose we perturb the particles away from their common circular 
orbit, but keep them in the same orbital plane. (Moving them out of the 
plane just introduces linearly independent vertical oscillation modes that 
do nothing of interest for us here.) If there were no other forces in the 
system, the particles would just execute epicyclic oscillations at frequency 
k about their original circular orbit. To make things more interesting, let 
us tether the particles together with a spring of spring constant K and zero 
equilibrium length, so that the spring always exerts a tension which in and 
of itself tries to pull the particles back together. Viewed in a Cartesian 
frame rotating at the original circular orbit angular velocity, with origin at 
the unperturbed particle position, s-axis in the radial direction and y-axis 
in the azimuthal direction, the linearized perturbed equations of motion of 
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the particles are easy to write down: 


-W = Xl ^- R -&R*' + m { * a ~ X,) ' 
d 2 2/i on d ^ K 

= “ 2fi_ rr + — ~ Vl )> 

at z at m 


(4.6) 


(4.7) 


and two similar equations for the accelerations of the second particle which 
are identical to these except with the indices 1 and 2 flipped. The first terms 
on the right hand sides of these equations represent the Coriolis acceleration, 
the last terms represent the spring accelerations, and the middle term in 
equation ([h(]) is the result of the combined gravitational and centrifugal 
acceleration. 

It is easy to find the normal modes of these coupled oscillator equations. 
Adding corresponding equations together, we get 

d 2 d clfl 2 

^2 ( x i + x *) = + y 2 ) _ r ~or^ Xi+ x 2 )’ ( 4 - 8 ) 

^ 2 ( 2/1 + 2 / 2 ) = + x 2 ). (4.9) 

Assuming a time dependence oc exp(— iujt), we immediately find that the 
oscillation frequency u> is given by 


w 2 = R- 


.dfl 2 

dR 


■ 4fl 2 = 


(4.10) 


which is hardly surprising. The particles are moving together as one, with 
an unstretched spring, executing epicyclic oscillations. (There is also a zero 
frequency mode corresponding to a constant azimuthal displacement of the 
particles.) The other normal modes are far more interesting. Subtracting 
corresponding equations, we get 


T 72 (24 - x 2 ) = 20-i(t/i - y 2 ) - R^-(x 1 - x 2 ) - — (x 1 - x 2 ), (4.11) 
d t z at dR m 


d 2 , d . . 2K 

T 72 (2/i - 2/2 = -2S2 — (xi - x 2 ) -(2/1 - 2/2)- 

d t z at m 

Hence the oscillation frequencies for these normal modes satisfy 

0 = w 4 - 


/ 4AT A 

, 2 2 K , 
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V m ) 
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V m 

d R ) 


(4.12) 


(4.13) 


This immediately implies that one of the solutions of this quadratic equation 
for lo 2 will be negative, implying an unstable/damped pair of frequencies, 
if and only if 
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All accretion disks around black holes have angular velocities that decrease 
outward. Hence a sufficiently weak spring always catalyzes an instability in 
the particle motion. Only if the spring is strong does it snap the particles 
back to their equilibrium position. 

In it simplest form, the MRI is exactly like this mechanical problem. 
Consider an equilibrium, differentially rotating fluid flow in which there is 
a uniform vertical component of magnetic field, as well as an azimuthal 
(toroidal) component which, while nonuniform, is too weak to exert signif¬ 
icant stresses on the equilibrium flow. Like both r- and g-modes, the MRI 
involves essentially incompressible motions, so let us assume for now that 
the flow can be treated as completely incompressible (V • v = 0). Also, to 
keep things simple, let us neglect buoyancy forces for now by assuming that 
the equilibrium has uniform entropy. Now consider local (i.e. wavelengths 
much less than equilibrium length scales) perturbations with wavenumbers 
purely in the vertical direction. The incompressibility condition then imme¬ 
diately implies that Sv z = 0, i.e. the fluid motions are entirely horizontal. 
Similarly, the condition V • B = 0 implies that there is no perturbation in 
the vertical held component, SB, = 0. 

The perturbed flux freezing equations are 


dSB R _ D dSv R 
dt z dz 


(4.15) 


and 


dSBtfi 

dt 


= R^SB r + B z 
d R 


dSv ^ 
dz 


(4.16) 


These equations describe the dynamics of the held lines (the spring!) them¬ 
selves. Vertical gradients of radial and azimuthal velocity stretch the equi¬ 
librium vertical held out into radial and azimuthal held, respectively. In 
addition, the background differential rotation shears out the perturbed ra¬ 
dial held into azimuthal held. 

The perturbed fluid momentum equations are 


dSv R 

dt 


— 2 QSv'f, 


B z dSB R 
Anp dz 


(4.17) 


and 


dSv ^ k 2 _ B z dSB$ 

dt 212 R Anxp dz 

Q dSP B 0 dSB 4 

dz A-np dz 


(4.18) 

(4.19) 


The radial and azimuthal momentum equations (4.17) and (4.IS) are very 
similar to the spring equations ( |4.6| ) and (4/7), the differences being due to 
the fact that we are now using Eulerian rather than Lagrangian variables. In 
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terms of the Lagrangian displacement £, the Eulerian velocity perturbations 
are 


<5v = 


. x-Tt; i: X7 

—+v-V£-£-Vv = —R + 


dt 


,an\ 


-Z« R lm )*■ < 4 ' 20 > 


Note that vertical gradients in the radial and azimuthal held are producing 
tension forces in equations (4.17) and ( 4. 18| ) very analogous to the spring 
forces considered earlier. The vertical momentum equation ( 4. 19| ) is rather 
interesting. Vertical gradients in the azimuthal held create vertical mag¬ 
netic pressure gradients which, for an incompressible huid, are immediately 
balanced by vertical huid pressure gradients. This equation just determines 
the huid pressure, and is otherwise decoupled from the other equations. 


Note that in the linear system the vertical pressure balance equation (4.19) 
becomes trivial if the equilibrium azimuthal held B$ vanishes, and in fact 
this is the only place where B$ appears. (We will see in section 5 that 
radiative diffusion can affect the MRI through this vertical pressure balance 
equation, however.) 

Assuming a space-tim e dep e ndenc e of th e pert urbat ions oc exp[j(&z — 
cot)], the four equations ( 4.15 ), ( 4.16 ), ( 4.17 ) and ( 4. IS ) can be solved to 
give the dispersion relation 

0 = u) 4 — (2kvA + k 2 ) w 2 + kv a (^kvA + > (4-21) 


where va z = B z /(Aicp) 1 / 2 is the A lfven speed correspo nding to the vertical 
held component. Equation (4.21) is identical to eq. ( 4.13| ), if we simply 
replace the natural frequency 2 K/m of the spring with ki’A z! 

It is easy to show from eq. (4.21) that the maximum instability growth 
rate is 


l^lmax — ty B 


an 


d R 


(4.22) 


and occurs for wavenumbers satisfying 


{kVAz)\u\max — 


16S7 2 ’ 


(4.23) 


Note that this is all in agreement with our statement above that r-mode 
physics will be substantially modihed when the natural frequency of our 
magnetic “spring” kv a is of order k. 

The analysis we have done so far recovers the simplest form of the MRI. 
It turns out that this mode is also an exact solution of the nonlinear local 
MHD equations, although it is itself vulnerable to further instabilities [ fl7| . 
It is often called the “channel solution”, as it is characterized by rapid 
counterstreaming motions in horizontal planes p^,|53|. 
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One can of course generalize the analysis to include general axisymmetric 
wavenumbers k = fc^R+fc-z and arbitrary equilibrium field directions. One 
can also add buoyancy by introducing entropy gradients. Most relevant to 
geometrically thin disks are vertical entropy gradients. In this case the MRI 
dispersion relation can be written as |l(| 

0 = N 2 ^ Co 2 - 4H 2 (k ■ v A ) 2 , (4.24) 

where Co 2 = w 2 — (k ■ v A ) 2 and v A = B/(47rp) 1 / 2 is the vector Alfven speed 
corresponding to the equilibrium field B. Note that magnetic fields are 
therefore directly connected to both r- and g-modes. The linear analysis of 
the MRI can also be extended to localized, nonaxisymmetric perturbations, 
which are again subject to unstable growth |ll|. A fully general relativistic 
linear analysis of the MRI has recently been done by Araya-Gochez 

The general local stability criterion in ideal MHD for rotating, axi sym ¬ 
metr ic fl ows turns out to be the same as the Hpiland criterion [eqs. (|4.l| ) 
and ( |4.2| )1, but with angular momentum gradients replaced by gradients in 
the angular velocity Vt |Hj : 

act 2 

IV 2 + R-irpr > o (4.25) 

oK 

and 

(VP x R) • (VH 2 x VS) > 0. (4.26) 

Magnetic fields appear to facilitate the existence of instabilities which are 
more directly tuned into the thermodynamic sources of free energy in the 
flow. Indeed, Balbus |H|] has shown that in situations where heat conduc¬ 
tion exists and is restricted to flow along magnetic field lines, the entropy 
gradients in the MHD Hpiland criterion are replaced by temperature gra¬ 
dients! This situation arises in cases where charged particles conduct the 
heat and the collisional mean free path exceeds the Larmor radius of the 
particles. In a non-rotating, hydrostatic system, the Schwarzschild stability 
criterion that the entropy increase upward is replaced by the criterion that 
the temperature increase upward! While interesting, it should perhaps be 
remembered that any temperature or angular velocity gradient constitutes 
a thermodynamic source of free energy, but magnetic fields are destabilizing 
for only certain directions of these gradients. 

The MRI relies on weak magnetic fields. If the field is too strong at the 
wavelength considered, then magnetic tension will overcome the effects of 
magnetic torquing and the flow will stabilize. The stronger the field, the 
longer the wavelength required for instability. This immediately suggests 
that laminar accretion disk flow models will be unstable provided an un¬ 
stable MRI wavelength ~ va/Q can fit inside the vertical thickness of the 
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disk. This corresponds to initial field energy densities that are less than 
the thermal pressure. That this does indeed provide an upper limit to the 
field strength in order for the MRI to exist is confirmed by global linear 
analyses fi(||45| . 

Since its discovery, many numerical simulations have been done by a 
number of groups that confirm that the nonlinear development of the MRI 
in an initially weakly magnetized medium leads to sustained MHD turbu¬ 
lence in the accretion flow. Although the simplest MRI modes are axisym- 
rnetric, it is essential that these simulations be fully three dimensional. Lo¬ 
cal axisymmetric simulations with an initial net poloidal field are inevitably 
dominated by the channel solution, which tends to break up in three dimen¬ 
sional simulations. Axisymmetric simulations with no initial net poloidal 
field only develop a transient phase of turbulence that eventually decays 
away. This is a consequence of Cowling’s anti-dynamo theorem, which in 
one version states that magnetic fields cannot be sustained by fluid motions 
under conditions of axisymmetry [ |S3| . 

Three dimensional nonlinear simulations of MRI turbulence have been 
done in a number of geometries, the simplest being the “shearing box” 
|n],|5^]. Here one zooms in on a tiny region of the disk and explores the 
local nonlinear behavior of the MRI, neglecting background gradients in the 
initial equilibrium except in the radial direction where one has to take into 
account the all important differential rotation. Periodic boundary condi¬ 
tions are adopted in the vertical and azimuthal directions, and also in the 
radial direction in coordinates that shear with the background differential 
rotation. (An excellent discussion of the shearing box can be found in [f57j|.) 
These shearing box simulations show clearly that the nonlinear development 
of the MRI is one of fully three dimensional, anisotropic MHD turbulence 
that exhibits strongly correlated fluctuations in azimuthal and radial com¬ 
ponents of velocity and magnetic field. Angular momentum is transported 
outward by a sum of Reynolds (fluid) and Maxwell (magnetic) stresses, i.e. 

tr^, = pv R 6v ^ - BR ^ 4 ' i (4-27) 

where Sv^ is the azimuthal velocity minus the mean background orbital 
velocity. The Maxwell stress generally dominates the Reynolds stress by 
factors of 3 to 4 |57|], and their values averaged over the computational 
volume both exhibit substantial variability on the orbital time scale. 

Unfortunately, shearing box simulations cannot tell us the magnitude of 
the anomalous stress, even in an averaged sense, appropriate for accretion 
disks in nature. The neglect of vertical stratification and buoyancy prevents 
the gas pressure from having much dynamical significance in the turbulent 
state. In fact, shearing box simulations carried out with an adiabatic equa¬ 
tion of state have monotonically growing pressure generated by heating in 
the turbulence. However, this increasing pressure does not affect the average 
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value of ^7|, perhaps because it acts simply to enforce incompressibility 
in the turbulence to a greater and greater degree. Instead of the pressure, 
the anomalous stress is more closely related to the overall magnetic energy 
density in the turbulent state, but this turns out to depend on the size of 
the computational box, the initial magnetic held strength if there is a net 
mean held in any direction in the box, the resistivity (numerical or physical) 
and the size of the artihcial viscosity used in the simulations |^,|5^]. The 
latter two quantities determine in particular the rate at which magnetic held 
“reconnects” in the simulations. The initial held topology also affects the 
level of turbulent transport. Simulations with no initial net poloidal held 
result in a stress that is roughly an order of magnitude below that obtained 
when there is a net initial poloidal held J57|,|5^|. 


Simulations have also been done in vertically stratified shearing boxes 


disk, fully taking into account the effects of vertical buoyancy and a hnite 
scale height of the disk set by vertical hydrostatic equilibrium. This latter 
fact is important, as it connects the local pressure to the vertical size of 
the disk. In simulations that start out with a weak magnetic held with 
no net poloidal component, three dimensional MHD turbulence is again the 
generic result, at least within the disk interior. The time and space averaged 
anomalous stress generated by this turbulence turns out to be proportional 
to the pressure at the midplane of the disk, in excellent agreement with the 
Shakura-Sunyaev prescription (HI)- The reason for this appears to be that 
the midplane pressure sets the vertical scale height. The thermodynamics 
of these simulations is relatively crude: in simulations run with an isother¬ 
mal equation of state, so the gas cannot heat or cool, the turbulent stress 
saturates to a constant average value. In simulations run with an adia¬ 
batic equation of state, so the gas continually heats as a result of turbulent 
dissipation, the scale height of the disk slowly increases in response to the 
heating, and the average turbulent stress also increases in lock step with 
the increasing midplane pressure |110[ . In agreement with the shearing box 
simulations, the stress is dominated by magnetic fields rather than velocity 
fields. 


11C, SOT . Such simulations zoom in on a small vertical slice through the 


The actual value of the proportionality constant a , as defined by 
divided by the midplane pressure, is typically ~ 10~ 2 , although different 
simulations can produce somewhat larger or smaller values 


11C, 


It is 

still not clear that this value is to be fully trusted, and it should certainly not 
be assumed to be a dimensionless constant of nature! In all the simulations 
that have been done so far, it appears that vertical buoyancy of the magnetic 
held does not play a dominant role in setting the saturation level of the 
turbulence, although this conclusion might be altered once the effects of 
radiative diffusion are incorporated (see section 5 below). Instead, local 
balance between magnetic held amplification by the MRI and dissipation 
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appears to be the most important factor in setting the overall stress level, 
just as in the shearing box simulations. Unfortunately, the dissipation is 
often largely numerical in nature, and the actual value of a produced by 
the simulation can depend on the grid resolution. (For example, doubling 
the grid resolution in one of the most recent simulations fiof increased the 
value of a by 1.5.) 


One dramatic effect that appears to be a generic feature of all the strat¬ 
ified shearing box simulations is the formation of a strongly magnetized 
corona above the disk [ 32, 110 ,8C |. This has been explored most extensively 
by Miller and Stone [flOf who considered an isothermal disk (resulting in a 
vertically Gaussian profile in pressure and density) within a computational 
domain that covered plus and minus five scale heights H z around the disk 
midplane. Figure |d| depicts the magnetic field lines in the turbulent state 
in one of their simulations. Within two scale heights of the midplane, MRI 
turbulence dominates and the field has a highly chaotic structure. Above 
two scale heights, however, the field becomes much more coherent as the 
field energy density dominates the gas pressure in the magnetized corona. 
Approximately one quarter of the energy generated within two scale heights 
of the midplane is transferred to the corona in this particular simulation, 
mostly in the form of Poynting flux associated with buoyant magnetic field 
lines. MRI turbulence therefore appears capable of self-consistently generat¬ 
ing an energetically important corona, in agreement with the observed fact 
that the hard X-rays in black hole sources can carry a substantial fraction 
of the total accretion power. 


Stratified shearing box simulations have also been done with an initial 
net poloidal field, but the results of the most recent such simulations produce 
a catastrophic change to a magnetically dominated disk structure that is 
inconsistent with the shearing box assumptions [Q. Large scale, ordered 
magnetic fields that might help generate an MHD outflow or jet also cannot 
be formed within the confines of a narrow shearing box. Global simulations 
of the entire disk are therefore of great interest, and some have now been 
completed. Examples can be found in |d],|7(||, 109 . Such simulations are 
quite challenging, and have their own problems when it comes to making 
inferences about how real flows in nature work. Because they attempt to 
simulate the entire disk, they are necessarily limited in spatial resolution 
and total duration in time. They also tend to suffer from transient effects 
because of the presence of an outer radial boundary condition or an imposed 
finite radial extent of the disk itself. Nevertheless, these simulations have 
generated some rather interesting results. 


First, time and space-averaged values of the Shakura-Sunyaev alpha 
stress parameter range from ~ 0.01 to ~ 0.2, consistent with results ob¬ 
tained from the local simulations. The reader should once again bear in 
mind all the caveats, however. Second, strongly magnetized coronal re- 
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Fig. 6. Magnetic field line structure generated by the MRI in a small section of 
an accretion disk, from a simulation by Kristen Miller. The figure shows a slice in 
the Z(vertical) - Y (azimuthal) plane of the disk (X is the radial coordinate). The 
color of the lines represents the strength of the field locally, normalized by the 
initial pressure at the disk midplane. MRI turbulence near the midplane of the 
disk (Z< 2 H z ) results in tangled small-scale field structures. In the lower density 
(Z> 2 H z ) regions, the field lines become smoother and show large scale coherence. 
(This figure was kindly generated by Kristen Miller. The field lines in this plot 
were generated by starting at points in a plane parallel to the X-Z plane and then 
integrating the magnetic field line vectors. This integration eventually becomes 
dominated by numerical uncertainties, particularly in the inner regions of the disk 
where the field is so chaotic, at which point the computation of the field line stops 
for this particular plot. Hence the appearance of a nonzero divergence in the field 
is an artifact of the way the field lines were computed: V ■ B — 0 is satisfied in 
the simulations.) 
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gions are generally formed outside the disk. Third, outflows are generated, 
although it is not yet clear whether a disk that starts out with a weak mag¬ 
netic field and develops MRI turbulence can generate a globally ordered 
field in the corona that can collimate and accelerate jets. There are re¬ 
cent intriguing results where rotating conical outflows that are externally 
confined by magnetic pressure are produced near the rotation axis fil|. [pif , 
but whether such outflows can produce the apparently self-collimating jets 
we observe on large scales in nature is far from clear. Fourth, significant 
stresses are exerted on the flow near the ISCO, a point that I mentioned in 
section 2 and to which I will return shortly. 

Before doing that, I would like to first address the issue as to whether and 
how far all this MRI physics can be used to justify the standard accretion 
disk equations that we wrote down in section 2. This has been carefully 
examined in the literature by Balbus and collaborators [|l5],[nj in the context 
of geometrically thin accretion disks, and I have made an attempt here to 
generalize their approach to include the effects of radial pressure support 
and advection. Our objective here is to see how far we can sweep all the 
MHD physics into a simple anomalous stress prescription. 

Using cylindrical polar coordinates (), the standard MHD equa¬ 
tions can be manipulated to express local conservation laws in mass, 

§ + V-(pv)=0, (4.28) 

radial momentum, 


d_ 

dt 


(, PVr)+V■ 


(pVRV - 



angular momentum, 
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dRj 4t tR’ 
(4.29) 
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and energy, 


d_ 

dt 


\pv 2 +P^ + pU+^) + 


V- 
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-pv 2 + pd> + pU + P 


v + 


— B x (v x B) 
47T 


+ F 


0(4.31) 


Here B p = BrR + B z z is the poloidal piece of the magnetic field and F is 
the radiative heat flux. 

There is nonaxisymmetric, time-dependent turbulence that is very com¬ 
plicated, but presumably there is some way of averaging over this turbulence 
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to get equations describing macroscopic, mean flow properties. It turns out 
that it is possible to (nea rly) r ecover the vertically integrated radial disk 
structure equations Q-iH) by such averaging, provided we make the 
following assumptions. First, the flow velocity is dominated by rotation, on 
top of which there are velocity fluctuations with root mean square ampli¬ 
tudes that are much larger than the mean radial infall velocity: 

v = Rflcf) + u with ^>< u 2 > 1 ^ 2 ^> | < ur > |. (4.32) 


Also, the turbulence is MHD in character and highly subsonic, with 

(j— ) ~< u 2 >«< c 2 >~ 1^. (4.33) 

\ 4t Tp / s < p> 


In addition, we assume that there is negligible mean magnetic field, i.e. 

| < B > | « B 2 > 1 / 2 . Finally, we assume that negligible fluxes of mass, 
angular momentum, and linear momentum leave the upper or lower disk 
surfaces, and the only flux of energy leaving these surfaces is radiative in 
character, i.e. contained in the radiation flux F. We already know that 
this last assumption is dangerous, as simulations clearly show substantial 
energy transfer into the disk corona. Provided it is only energy and not 
mass, angular momentum or linear momentum, we will be able to get away 
with this by calling F an energy flux as I did in section 2. However, this 
immediately precludes the possibility of developing large scale outflows that 
carry away mass and angular momentum from the disk. 

Now, following & define a vertically integrated, density weighted, spa¬ 
tially averaged value of a quantity X at radius R in the disk by 


< X >= 


Mpx) 

s 


i 

2ttART, 



r R+AR/2 


dR 


lR-AR/2 



(4.34) 


where A is a vertical integration and averaging operator, A R ~ H <C R 
is a small radial scale length over which we are averaging (assuming the 
turbulence causes variations on size scales smaller than a measure of the disk 
half-thickness H ), and X = Ap is th e surface density. Applying the operator 
A to the continuity equation ( 4.28 ), and using the boundary conditions of 
vanishing mass flux on the disk surface and azimuthal periodicity, we obtain 


I r) 

ar + 5^ E< "*» = 0 - 


(4.35) 


where the radial derivative is now defined by differencing over the range 
(R — AR/2,R + AR/2 ). Now, by a stationary flow I mean one in which 
the average quantities I have constructed are independent of time. Then 
i?X < ur > is a constant, which I may write as 


M = —27ri?X < ur > . 


(4.36) 
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Using the same procedure on the rest of the conservation equations ( 4.29 ) 
(4.31), I obtain 


ri P 


■ de d 
M d R ~ d R 


and 

M 

where 
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d R dl? \p J 
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I have been forced to define other averages here, in particular 

< VrU > 


U = 


<u R > 


and 
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P = < Pvr/p > 
P <u R > 


(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 


The last two equations may not be rigorously consistent, but this is the 
usual consequence of attempting to describe a two-dimensional average flow 
with one-dimensional, vertically integrated equations. 

If I wri te E = 2 p H, then equations ( ff.36| )-( |4.39| ) become identical to 
equations ( |2.7| )- (|2. 1C ) that we discussed in section 2, with one important 
exception. Note that we recover the true meaning of the anomalous stress 
in equation ( |4.4C| ) : the correlated radial and azimuthal velocity and mag¬ 
netic field fluctuations that we saw in equation (4.27)! We have (almost) 
succeeded in hiding the underlying MHD physics in the average flow by 
placing it entirely inside the anomalous stress. 

The only problem is that I w as unable to r ecove r the left hand side of the 
radial momentum equation (2.8) in equation (4.37), i.e. the term pvdv/dR. 
The reason is that I chose to neglect the magnetic terms in the original 
exact equation ( 4.29 ). I am justified in doing this, and still retaining the 
radial pressure gradient term, because the turbulence is subsonic. The same 
reasoning, however, forces me to neglect the poloidal velocity terms in this 
equation. If I want to retain them, then I must also retain the magnetic 
terms. 
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In hydrodynamic models of accretion disks with radial advection, the 
pvdv/dR term only becomes important near the sonic point, but it is pre¬ 
cisely here that the assumptions underlying the averaging process giving 
rise to equations (4.36)-(4.39) fail utterly. The poloidal velocities are tran¬ 
sonic and the character of the underlying turbulence must change. Mag¬ 
netic fields are likely to be important here as they are stretched out and 
sheared by the large inflow velocities, and we cannot neglect them. It is this 
physics that is the source of inner torques exerted on the disk across the 
ISCO |^,Q. Rather than localized turbulence over which we can average, 
we have torques that are generated by more coherent magnetic field lines 
formed by flux freezing in the rapidly infalling fluid. Averages cannot be 
performed here because the time and space dependence of the “turbulence” 
is comparable to that of the macroscopic flow itself. 


Numerical simulations that probe the nature of these non-MRI inner 
torques have been done by a number of authors |||, [l(J. The most 
recent such simulation, which has the highest resolution, confirms the ex¬ 
istence of significant stresses in the plunging region inside the ISCO, with 
matter having ten percent more binding energy and ten percent less angular 
momentum ||6(j|| . The Shakura-Sunyaev a parameter rises to huge values in 
the plunging region, but this is a mis-application of that prescription be¬ 
cause it is due merely to the drop in pressure in this region of the flow. 
The flow is highly time-dependent and cannot be characterized by station¬ 
ary models flowing through critical points of ordinary differential equations. 
As always, there are still caveats in these simulations and further work is 
needed. In particular, all the simulations that have been done so far lack a 
full general relativistic treatment. Performing a simulation in a Kerr space- 
time in which the black hole can exert magnetic torques on the accretion 
disk, thereby supplementing accretion power with spin power, is likely to 
be a profitable research direction! 


I have focussed in this section on angular momentum transport in black 
hole accretion flows, but the existence of turbulence also begs the question 
as to whether or not there could be substantial anomalous heat transport 
as well. Unfortunately, none of the simulations that I discussed here can 
address this question, as they do not include radiation transport. In fact, 
the thermodynamics of the fluid is generally treated very crudely by simply 
adopting an adiabatic (P oc p 5 / 3 ) or isothermal (P oc p) equation of state. 
Dissipation in the turbulence is not properly accounted for, and energy is 
generally not conserved. It turns out that incorporating radiation transport 
directly into the physics of the MHD turbulence leads to very interesting 
results, and this will be the topic of the next and final lecture. 
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5 The Role of Radiation Magnetohydrodynamics 

Shakura-Sunyaev based models of standard accretion disks around black 
holes are always radiation pressure dominated if the disk luminosity is any¬ 
where near Eddington. If we take r ^ = aP, with P being the total (gas 
plus radiation) pressure, and assume that vertical heat transport proceeds 
through radiative diffusion, then the ratio of radiation to gas pressure is 
roughly 


Trad 


— 9 <w,,, 


where / is the fraction of local accretion power that is not dissipated in 


the disk interior 11141. As I discussed previously, there has always been 
uncertainty as to whether the anomalous stress scales more with radiation 
pressure or gas pressure, or indeed something else, and this has led to ambi¬ 
guity about the thermal and “viscous” stability of the inner parts of black 
hole accretion flows. Even if we take = aP gas , which is thermally and 
“viscously” stable, the ratio of radiation to gas pressure is still large: 


Trad 

pT 


gas 


io 4 a‘/V s/5 (i - /) 9/5 (j^) 


1/5 


(A-f 7* 

V bEdd / 


- 21/10 
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(5-2) 

Note that radiation pressure is most important for models of active galactic 
nuclei and quasars that contain supermassive black holes, but it can still 
be important for stellar mass black hole X-ray binary models. As noted in 
section 3 above, allowing much of the locally dissipated power to vertically 
escape from the disk interior in non-radiative fashion (/ —> 1) reduces the 
importance of radiation pressure [114]. 

In addition to the issue of secular stability, there are also other uncer¬ 
tainties in this inner portion of the flow. Radiation pressure support against 
the tidal field of the black hole implies that 


n e a T p ^ Q^p Z . ( 5 . 3 ) 

c 

The electron number density n e ex p, so that the mechanical equilibrium 
condition ( |5.3| ) gives F (x z, on dF/dz equals a constant. On the other 
hand, it has long been assumed that the vertical distribution of turbulent 
dissipation per unit mass is probably constant in the disk interior, i.e. that 
the power dissipated per unit volume is proportional to the local density. 
Radiative transport of heat then implies that 


V F ~ 


dF 

d7 


oc p. 


(5.4) 
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Hence mechanical and thermal equilibrium immediately lead to the conclu¬ 
sion that the density p is independent of height z. Such an equilibrium 
is violently unstable to vertical convection |2^]. For example, eliminating 
the density gradient in equation (4.4) leads to an imaginary Brunt-Vaisala 
frequency. More physically, consider an upward adiabatic perturbation of 
a fluid element in a constant density background. Because the pressure 
drops with height, the fluid element expands and becomes less dense than 
the surrounding constant density background. It is therefore buoyant and 
unstable. 

We are therefore lead to the conclusion that heat is transported out of 
the disk convectively, not radiatively. However, this whole argument is inti¬ 
mately wrapped up in the physics of the MRI. As noted in the last section, 
the MRI modes are tightly coupled to the gravity modes that are responsible 
for hydrodynamic convection. How does convective heat transport work in 
the presence of MRI turbulence? Moreover, the convective instability arises 
because of our assumptions about the mass distribution of dissipation, and 
this too is determined by the MRI turbulence, so how does this really work 
in a vertically stratified, radiation pressure dominated medium? 

It turns out that radiation MHD profoundly affects the MRI, the Parker 
instability, and also introduces new types of dynamical instability in the 
flow. Radiation MHD is just like standard MHD, except that the photon 
“gas” is treated as a separate fluid that couples to the plasma through 
absorption, emission, and scattering. I summarize the basic Newtonian 
equations of radiation MHD that are useful in studying the physics of the 
radiation pressure dominated regions of high luminosity black hole accretion 
flows in the appendix. The reason that we have to treat the photons in a 
more careful way than just writing the total pressure P to t as the sum P gas + 
P ra d is that, in addition to being effectively thin (r e ff < 1) in some cases, 
the innermost regions of black hole accretion flows do not have tremendous 
total optical depth either. For = aP ra d, the vertical Thomson depth of 
the disk midplane is 


tt ~ a Vl - /) 


- IT 2 
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(5.5) 
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On the other hand, tr^ = aP gas produces a much denser, more optically 
thick disk, 


tt ~ 10 4 a 4 / 5 ?7 3 / 5 (l — /) 1//o 
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Even in this case, photons, which are providing the dominant pressure sup¬ 
port, are highly diffusive, especially on small length scales. 

Even in a perfectly electrically conducting fluid, significant radiation 
pressure with radiative diffusion alters MHD in two important ways. First, 
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the compressibility of the fluid is greatly enhanced on small scales beyond 
what would naively be expected based on a very high radiation sound speed. 
Even highly subsonic motions [speeds <C (4aT 4 /'dp) 1 /' 2 ] will be compressible 
if they occur on a scale small enough for photons to diffuse during the 
motion. Second, temperature fluctuations in the gas tend to be smoothed 
out by the radiation field. These effects have important consequences for the 
MRI and the Parker instability, and they also produce wholly new classes 
of instabilities. The nonlinear, turbulent state of the radiation pressure 
dominated region of an accretion disk is therefore likely to be very different 
from that envisaged by standard accretion disk theory, and the sorts of 
modeling we have been discussing up to now may be doing a very poor job 
of describing this, energetically most important, region of the flow. 

Radiation MHD effects on the Parker instability were considered quite 
early in the development of accretion disk theory J98| . Consider an isolated, 
straight, horizontal tube of magnetic flux immersed in the vertically strat¬ 
ified, radiation pressure dominated environment of a black hole accretion 
disk. If the tube has small diameter, radiative diffusion will tend to smooth 
out temperature (and therefore radiation pressure) differences between the 
interior and exterior of the tube. The excess magnetic pressure inside the 
tube must therefore be balanced largely by a deficit of gas pressure caused 
by a lower interior density. The flux tube therefore has a tendency to be 
buoyant. The standard Parker instability in pure gas disks implies that 
flux tubes will be buoyant if B 2 /(8 tt) > P , and this then suggests that 
tubes with diameters less than the radiation diffusion length over an un¬ 
stable mode growth time will be unstable if B 2 /(8n) > P gas . That this is 
in fact true is supported by numerical simulations of flux tube dynamics 
in radiation pressure dominated environments [ p9[ . Specifically, if the flux 
tubes are strong enough to self-consistently maintain an anomalous stress 
tri/) = a-Ptot) they are extremely buoyant. If instead the flux tubes are 
only as strong to maintain the lower stress tr^ = aP gas , then they can 
be retained in the medium for much longer time scales. This is at least 
partly due to the fact that the ambient medium is then much denser and 
more optically thick, so that photon diffusion across the flux tubes is much 
slower. Unfortunately, the results of this study are not conclusive, as they 
fail to take into account the MRI (whose application to accretion disks was 
discovered two years later). The question remains as to whether the MRI 
can either generate the field sufficiently fast to counteract rapid buoyant 
escape, or directly modify the buoyant dynamics itself. 


Radiation MHD also modifies the behavior of the MRI. Recall that the 
MRI is nearly incompressible. Even in the nonlinear regime, B 2 /(8 tt) is 
substantially less than the thermal pressure in the disk interior and the 
turbulent motions are highly subsonic. However, radiative diffusion will 
modify this. From equations (4.22) and (4.23), the characteristic time scale 
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associated with the MRI is the orbital period and the characteristic length 
scale is ~ va/Q- The turbulence will therefore remain incompressible only 
if the orbital frequency times the diffusion time tdiff over a distance v\/Q 
is substantially greater than unity. Now, 




(RT P\ 

(<) 

TR(f,K T 

V C ) 

U 2 J 

cfl 


(5.7) 


where the last equality comes from equation (4.40) and the fact that MRI 
turbulence is dominated by magnetic rather than Reynolds stresses. On the 
other hand, if heat is transported vertically by radiative diffusion and the 
disk is in thermal equilibrium, then equations (0) and (]2.12j) imply 


T R^ T Jfad 1 

cH ~ E 2 fi 2 ~ 


(5.8) 


the last equality coming from vertical hydrostatic equilibrium. Hence in a 
radiation pressure dominated disk in which heat is transported radiatively 
and the turbulent stress is of order the magnetic stress, then flfdiff ~ 1 and 
the turbulence is therefore expected to be compressible [119]. Note that this 
argument is rather general, and in particular is completely independent of 
the Shakura-Sunyaev stress assumption. 

Compressibility affects the MRI in interesting ways. Consider again 
the channel solution that we discussed in section 4. We noted that the 
equilibrium azimuthal field does not affect the linear growth of the channel 
solution at all. In fact this is true for all the linear axisymmetric MRI 
modes, whose growth depends on the equilibrium held only through k • va 
from equation ( (4.24[ ). If there is some equilibrium B^ (and if there is not, 
there soon will be!), then the channel solution rapidly develops vertical 
magnetic pressure gradients. Normally these do not affe ct an ything as the 
fluid pressure is able to balance these through equation ( |4.19| ). However, if 
the radiation is diffusive on these length scales, then the radiation pressure 
is lost. If the azimuthal held energy density exceeds the true gas pressure, 
then not even the gas pressure can balance the vertical magnetic pressure 
gradients. As a result, vertical motions must be excited which take energy 
out of the growth of the MRI. This effect does not kill the MRI, but it can 
dramatically reduce its linear growth rate . If the azimuthal held energy 
density exceeds the true gas pressure and if B^ > B z , then the growth rate 
of the MRI is reduced below the orbital frequency by the ratio of the gas 
sound speed to the azimuthal Alfven speed. This suggests that growth of 
magnetic stresses above the gas pressure may be sluggish, but numerical 
simulations are really required to answer this question. 

Such simulations have now been done in non-stratihed shearing boxes, 
both in two fllg|| and three fll9|| dimensions. Three dimensional simulations 
that have a net vertical magnetic flux through the box develop magnetic 
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stresses that are independent of radiation diffusion effects - while linear 
growth of the channel solution may become sluggish, it still grows! On the 
other hand, it appears that radiation diffusion does diminish the turbulent 
magnetic stresses that develop when there is no net vertical magnetic flux 
through the box. Unfortunately however, for the same reasons that we 
discussed in section 4, shearing box simulations cannot provide definitive 
estimates of the level of stress that will exist in real flows, and vertically 
stratified and global simulations including radiative diffusion will need to 
be done. 


Perhaps the most interesting result that has emerged from the recent 
shearing box simulations is that radiative diffusion does indeed make the 
nonlinear turbulence highly compressible, and large density fluctuations (by 
more than a factor 20 1191!) form and reform in a highly time-dependent 
manner. One important consequence of this fact is that some of the tur¬ 
bulent energy is directly dissipated into photon energy on scales that are 
completely resolved by the simulation. This occurs because photons diffuse 
out of compressed regions and into rarefied regions, damping compressive 
motions in the fluid Q. Hence not only are we now getting a handle on 
the behavior of the anomalous stress with respect to angular momentum 
transport in the radiation pressure dominated region, we are also learning 
about how the mechanical energy is dissipated into heat. Ironically, this is 
easier to do here than when gas pressure dominates! 


There are also entirely new classes of instabilities that are expected to 
be present due to radiation pressure and radiative diffusion when the ver¬ 
tical gravity is included. One of these is the “photon bubble” instability, 
first examined by Arons in the context of accreting X-ray pulsars |HJ] . Here 
density fluctuations in the presence of a strong magnetic field are driven 
buoyantly unstable by photons diffusing into and heating (out of and cool¬ 
ing) underdense (overdense) regions. A local MHD instability that may 
have similar physics has also been found to exist in the radiation pressure 
dominated region of black hole accretion disks . Fast and slow magne- 
tosonic waves at short wavelengths where radiative diffusion causes the loss 
of radiation pressure support of the wave can also be unstable to periodic 
driving by the equilibrium vertical forces (radiation pressure gradients, gas 
pressure gradients, and gravity), at least for waves propagating in certain 
directions |j(|. The characteristic length scale for these instabilities is the 
gas pressure scale height (tiny in a radiation pressure dominated medium) 
and the growth rate is much larger than the orbital frequency, by a factor 
of order the square root of the radiation to gas pressure ratio. These insta¬ 
bilities can only exist in Thomson scattering dominated media if magnetic 
tension forces help support the wave, which always occurs for wave vec¬ 
tors that are neither parallel or perpendicular to the equilibrium magnetic 
field. (In some cases, the small absorption opacity can drive even purely 
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hydrodynamic acoustic waves unstable Q.) 

Simulations have yet to be done that track the nonlinear development 
of these instabilities in an accretion disk environment, but Begelman |2^| 
has constructed a one-dimensional, nonlinear periodic shock train solution 
that may describe the outcome of the magnetosonic wave instabilities. This 
solution produces a medium that is highly dynamic and porous to radiation 
diffusion, and can in fact support fluxes that are super-Eddington by factors 
of 10-100 without driving an outflow J24| (see also [ 1051)! 

Both the density fluctuations generated in the shearing box simulations 
of the MRI and these compressive wave instabilities strongly suggest that 
radiative diffusion generates a highly inhomogeneous, very time dependent, 
somewhat radiatively porous inner accretion flow. This speculation needs 
to be verified by further work, but if true, there are numerous observational 
implications, including possibly super-Eddington fluxes and modification of 
local spectra due to increased thermalization of the photons with the denser 
phases of the medium. The standard reasoning behind the putative exis¬ 
tence of thermal/“viscous” instabilities in this part of the flow will also be 
affected. Suppose for example that a perturbative increase in temperature 
increased the porosity of the medium so that photons could escape more 
easily. Such a process could stabilize the disk against a thermal runaway. 
It has even been suggested that the radiation pressure dominated region of 
the disk may break up into actual discrete dense clumps interpenetrated 
by a much hotter, more optically thin medium capable of producing hard 
X-rays by Comptonization fl6|| . Clearly, further numerical simulations are 
likely to produce very interesting results! 


6 Conclusions 

After these lectures, particularly the last one in section 5, a student may be 
feeling a strong urge to run screaming from the incredibly complex physics 
that must be understood before we can truly build realistic models of ac¬ 
creting black hole sources. Let me urge the reader that this was not my 
intent. The physics is indeed complex (and fascinating), and this essential 
fact is the reason why it has proved so hard to move beyond the Shakura- 
Sunyaev stress prescription introduced almost three decades ago. However, 
this is not the time to leave the field, but to enter it! For the first time 
we are beginning to investigate the basic physical principles that underly 
these flows, and we also have (or at least will have in the very near future) 
the simulation hardware that can provide definitive answers to the concrete 
questions we are now posing. The time is now ripe for genuine and definitive 
theoretical breakthroughs that will lead to predictive models that can ex¬ 
tract truly useful physical information from observations of accreting black 
hole sources. I encourage you to join in and take part in these exciting 
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developments. 


I would like to thank Shane Davis and Kristen Miller for preparing some of the figures in 
these lectures, as well as my collaborators (Eric Agol, Ivan Hubeny, Julian Krolik, and 
Aristotle Socrates) for their scientific insights over the years. I also acknowledge very 
useful exchanges with Ian George, John Hawley, Pierre-Yves Longaretti, Patrick Ogle, 
and Neal Turner. I am grateful to the organizers of the Les Houches summer school for 
the opportunity to participate in this excellent program! 


A The Equations of Radiation Magnetohydrodynamics 


For completeness, I su mmar ize here the basic (Newtonian) equations of radiation magne¬ 
tohydrodynamics (e.g. ]l 1 1 | ) that are relevant to black hole accretion disk applications. 
These are the mass continuity equation, 

g+V-(pv) = 0, (1.1) 

the gas momentum equation, 

p(^ + vVv ) = -v P + pg + T ( v xB) x B + ^ F , (i.2) 

the total internal energy equation, 

— + v jVj(u + E) + (u + E)X/jVj = -pVjVj - - V jFj , (1.3) 

the radiation energy equation, 


dE 4 
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the radiation momentum equation, 


_L ( dFi 

c 2 \ dt 


+ v-j V jFi + jVj ] — V jPij — 


k f p 


Ei, 


and the flux-freezing equation, 


— = Vx(vxB). 
dt y ’ 


(1.4) 


(1.5) 


( 1 . 6 ) 


Equations Jl.t - diJ ) have been written in Cartesian tensor notation, with summation 
implied over repeated indices. Although all the equations above are Newtonian equations, 
v/c effects have been included in the radiation terms, and all radiation quantities are 
defined in the local fluid rest frame. The quantities p, u , p, and T are the density, 
internal energy per unit volume, pressure, and temperature in the gas, respectively. They 
are related by equations of state, 

(1.7) 


7-1 


where 7 = 5/3 is the ratio of specific heats for an ionized gas, and 


P = 


pkT 


( 1 . 8 ) 


where k is Boltzmann’s constant and p is the mean molecular weight of the gas particles. 
For simplicity, I have neglected the effects of ionization and recombination in the gas. 
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The quantity v is the gas velocity and g is the local gravitational acceleration presumed 
to arise exclusively from the central object (the black hole). 

The radiation energy density E , the radiation flux F, and the radiation pressure 
tensor Pij are defined as frequency-integrated angular moments of the specific intensity 
in the local fluid rest frame: 


E = 


Fi = 


f 


and 


dvEv = — 

/ dis (p dQ.Iv( n), 
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The equations can be closed by integrating the radiative transfer equation directly to 
solve for the tensor Eddington factor fij , defined by 


Pij — fij F • 


( 1 . 12 ) 


The various opacity factors (all in units of cm 2 /g) are defined in terms of different fre¬ 
quency averages over the radiation field. The energy mean, Planck mean, and flux mean 
opacities are defined by 


k b = — / dvrfj{p,T)E„ 
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(1.13) 


(1.14) 


(1.15) 


respectively. Here xj, h (p, T) is the thermal absorption coefficient (in units of cm -1 ) at 
frequency is, n e is the electron number density, and crp is the Thomson cross-section. 
The Thomson opacity is 

- neCr T n 

Kp = -. (1.16) 

P 

Local thermodynamic equilibrium in the gas has been assumed, so that the thermal 
source function in the gas is given by the Planck function B V (T ) at the gas temperature. 

Due to the fact that it dominates the opacity in the inner portions of black hole 
accretion disks, I have incorporated the effects of electron scattering in the equations. 
This includes Compton scattering, which can dominate the thermal coupling between the 
gas and the radiation in some cases. I treated this in the frequency diffusion (Kompaneets) 


limit by approximating the radiation field as isotropic, 


The quantity is in the radiation energy equation (1.4) is defined as a frequency average 


n the local rest frame (see e.g. | 


over the local rest frame radiation field in the following way, 


is = 


— / disisE u 



E,c 3 \ 
87 xhis^ J 


(1.17) 


(The second term inside the parentheses in this equation represents the effects of stimu¬ 
lated scattering, which is necessary to include if one wants zero heat exchange between 
the gas and radiation when they reach exact thermal equilibrium.) 

For a superb introduction to radiation hydrodynamics (sans magnetic fields!), see 
the book by Mihalas &; Mihalas Jf9||. 
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